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We use a light-cone approach to solve the Schwinger-Dyson equation for the quark pro- 
pagator. We show how this method can be used to solve the equation beyond the space-like 
region, to which one is usually restricted with the Euclidean-space approach. We work in 
the Landau gauge, and use an infrared-enhanced model for the gluon propagator and in- 
clude instanton effects to get both confinement and vacuum condensates. With our models 
reasonable fits to known quantities are obtained, resulting in a light-cone quark propagator 
that can be used for hadronic physics at all momentum transfers. 



I. INTRODUCTION. 



For a microscopic QCD description of hadrons and hadronic matter one needs the fully dressed non- 
perturbative quark and gluon propagators, for which the Schwinger-Dyson formalism is a natural approach. 
A full study of QCD, however, requires investigation of hadronic properties at all momentum transfers. 
Since instant form of field theory is difficult to use for composite states at medium or high momentum, a 
light-cone representation is desirable In the present paper we develop a light-cone formulation of the 
Schwinger-Dyson equation for the quark propagator for use in developing hadronic light-cone Bethe-Salpeter 
amplitudes as well as providing new aspects of the quark propagator, which wc discuss below. 

It is well known that short distance processes occurring inside hadrons are well described by perturbative 
QCD (p-QCD) calculations. It is also well known that quarks never emerge as such from the reaction region, 
but instead they hadronize at a distance of about a fermi. But often these processes are not well described 
by p-QCD. The reason is that perturbative calculations in QCD have a rather limited range of validity: they 
are a good approximation only at very short distances where the expansion parameter (the effective QCD 
coupling constant) is small. But if QCD is the correct theory of the strong interactions (and that is the 
consensus today), then it should describe medium and long distance hadronic processes as well. We need 
then to understand how to compute the basic entities of the theory (its Green's functions) not just over very 
short distances, but over long distances too. 

The Schwinger-Dyson equations of a field theory embody all its dynamics p. 475]. They are the 
complete equations of motion for the Green's functions of the theory, and thus provide a natural way for 
studying the theory beyond the limited scope of perturbative expansions. Unfortunately they consist of an 
infinite tower of coupled integral equations relating full n-point functions to full [n + l)-point functions. 
Thus, the integral equation satisfied by one propagator may involve another propagator and a three-point 
vertex. The equation for this vertex may in turn include another propagator and a four-point vertex or 
scattering kernel, and so on. Some physically motivated truncation scheme becomes mandatory before the 
infinite tower can be brought to a manageable size. As it has been stressed in S the Ward identities of gauge 
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theories significantly ease this truncation, since they imply that two-point functions (propagators) uniquely 
determine the longitudinal part of three-point functions (vertices). 

The Schwinger-Dyson equation for the fermion and gauge boson propagators in QCD and QED have been 
studied using different approximations and models. For an excellent review see Q and references therein. 
A recurring topic in this area is the question of the analytic structure of the propagators. For the electron 
propagator one expects a singularity at = rripf^y^ . For a confined particle like the quark it is not so clear 
what one should expect, it depends on what the confining mechanism is. Coleman pp. 378-386] has 
shown that it is perfectly possible for a confined particle to have a pole in its propagator at a positive real 
(i.e., physical and time-like) value of p^. 

Fermion propagators with branch points at complex conjugate locations off the real axis of the variable 
p^, have been obtained since more than twenty years ago for QED ||], and more recently for QCD as well ||]. 
This has been studied in great detail by Maris In the case of QED, where we know that there must 
be a singularity at — rUpf^y^, this is thought to be an artefact of the approximations with no physical 
significance. For QCD, the absence of singularities on the real axis has been thought to be related to quark 
confinement. §avkli and Tabakin Q have shown that the imaginary part of the value of p^ where the 
singularity occurs may be interpreted as a decay width for a single quark state, and can be related to the 
hadronization distance. However, the fact that in two physical situations as different as those of QED and 
QCD the singularities seem to have the same origin [zj makes any physical interpretation difficult. The 
location of these singularities poses yet another difficulty: it invalidates the Wick rotation, and thus makes 
the theory as defined in Euclidean metric not equivalent to that defined in Minkowski metric. 

In this paper we study the Schwinger-Dyson equation for the quark propagator starting from its Minkowski- 
space formulation. We do make assumptions about the location of the singularities similar to those necessary 
to justify a Wick rotation. However, we arrive at a formulation that allows us to solve for the quark 
propagator in the space-like region (which would correspond to a solution in Euclidean space) or, in principle, 
to extend into the time-like region as far as needed. 

Our derivation of the Schwinger-Dyson equation in a light-cone representation is discussed in detail in 
section P The method is an extension of the method used in perturbation theory j^] to obtain an infinite 
momentum frame formulation, which in many respects is equivalent to the light-cone field theory formulation 
of perturbative diagrams. It is also analogous to the method for obtaining a light-cone representation of the 
Bethe-Salpeter equation from the standard four- dimensional instant-form, however, as explained in section 
the simple prescription for projecting onto the light cone cannot in general be used for the SD equation. 

In section [II we discuss our models. In sections IV and ^ we discuss some technical aspects of solving 
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the equations and give sample results in Minkowski space. In section VI we discuss the results of our model 



calculations and in section VIl give our conclusions 



II. THE QUARK SD EQUATION AND THE LIGHT-CONE. 



In this section we discuss our light-cone formulation the Schwinger-Dyson equation. Since the light- 
cone representation of the Bethe-Salpeter (BS) equation for bound systems is described extensively in the 
literature, we also briefiy review the BS equation to motivate the present work and explain its applicability. 
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FIG. 1. Diagrammatic representation of the Schwinger-Dyson equation for the quark propagator. 



The full dressed quark propagator satisfies the Schwinger-Dyson (SD) equation (see Fig|l|): 

s~\p) = s^\p) - Up) 

(1) 

Up) =^9' J ^T^oaS{qmq,p)D;lip- q), 

where rj^((7,p) and D'^jj{p—q) are the dressed quark-gluon vertex and dressed gluon propagator. Greek letters 
represent Lorentz indices and Latin letters stand for color indices. The color structure of all quantities in 
(|l|) is known: 

The bare vertex: T^^ = 7^Aa/2 

The dressed vertex: r^(g,p) = T''{q,p)\b/2 

The gluon propagator: D'^^^ik) = 5''^Df,^{k) 

so that the color indices can be contracted: AaAf,(5°'' = 16. The inverse bare propagator is S^^{p) =^ — rric , 
with the current quark mass rric zero in the chiral limit or estimated from PCAC. 

The solutions of the Schwinger-Dyson equation are of the form 

S-\p) = A{p')^-B{p^). (2) 

Physically, the most significant aspect of the solution for the two functions A{p'^) and B{p^) is the ratio 
M{p'^) = B{p'^)/A{p'^), which is interpreted as the effective mass of the dressed quark propagator. Since 
isolated quarks are confined, the interpretation of this mass is not straight-forward as for the electron SD 
equation. 

The light-cone formulation of quantum mechanics Q starts with the demonstration that in the light-cone 
representation one obtains light-cone Poincare generators, with an interaction-free Lorentz boost in one 
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direction. One can also show that the analogous result is true in a light-cone field theory. This enables 
one to study high momentum transfer processes, such as form factors in the region above 1 GeV, which 
is very difficult in the instant form. Although there was work on the SD equation three decades ago for 
perturbative QED for nonperturbative QCD the problem is quite different, as we show below. Before 
discussing the light-cone representation of the Schwinger-Dyson equation we review the well-known light- 
cone representation of the Bethe-Salpeter equation. Since in physical application of the BS amplitudes to 
form factors and transition amplitudes one needs the dressed quark propagator, the main motivation of our 
present work is to formulate a light-cone Schwinger-Dyson (LCSD) equation to be used with the light-cone 
Bethe-Salpeter (LCBS) formalism. 

There has been a great deal of work on bound state problems in light-cone representations based on the 
Bethe-Salpeter (BS) equation. For example, the BS equation for the pion in relative coordinates has the 
form 

= J dHjC{k,l)-^{l), (3) 

where /C is the kernel and \E' is the BS amplitude. One can obtain a light-cone representation [|lO| by inserting 
a delta function to include the light-cone on-shell condition and eliminate the l~ variable 

*(x,ki)= y[dy][d2lJ/C(x,ki;y,U)vI/(y,l^). (4) 

BS amplitudes from this type of equation have been used to study the transition from nonperturbative to 
perturbative regions We emphasize that for such calculations both on-shell and off-shell aspects of 

the SD functions A{p^, B{p^) are needed. Therefore the prescription of projecting onto the light cone as 
for the BS equation is not appropriate.. We now discuss our approach to the the SD equation for the quark 
propagator on the light cone. 

We solve equation (|l|) in a light-cone representation by using a method introduced originally for perturba- 
tion theory. In Ref j^S. Chang and S. Ma showed how the rules of light-cone perturbation theory (LCPT) ^ 
can be derived from the usual covariant Feynman rules simply by changing into light-cone variables: 

In simple calculations, like a one-loop self-energy diagram, an important feature of the new rules becomes 
apparent: the range of integration over the variable becomes finite 

+ 00 <.p+ 

dq+ dq+, 



where is the "plus" (longitudinal) component of the external momentum. As shown in that paper, this 
feature is related to the fact that in LCPT there are no diagrams with lines going backwards in time and no 



'^More accurately, their work refers to the Feynman rules in the infinite momentum frame, but the difference is of 
no bearing in our present discussion. 
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vacuum diagrams (with the exception of zero modes and instantaneous terms in fermion propagators). This 
is a crucial property of hght-cone field theory, since it simplifies tremendously the structure of the vacuum. 

In solving the Schwinger-Dyson equation we are faced with the non-perturbative self-energy diagram seen 
in Fig. |l|: 



I 



q^^q^^ . . . q^'^fQ {q^ + ze) fa {{v - qf + ^e) (5) 



Consider first the scalar case n = 0. Using the variables: 

a = q^ Iv^ ■, s' = q - q, s = p ■ p, q\ = - a pj 



the integral in becomes: 



^/q (s' + ie) fa [^a-^V{ql, s, s' , a, e)) , (6) 



2|a|(27r)« 

where all the integrals are from minus infinity to plus infinity, and 

Viq^, s, s', a, e) = q^ + a{l — a){—s) + (f — a)s' — iae 

is reminiscent of the similar quantity that appears after the well known Feynman trick of combining denom- 
inators. Particularly, if s' — tjiq — ie and /g(s) — ~ then after the q'j^— integration 



{a{l — a){—s) + (1 — a)mQ + amf — ie) 



a very familiar expression. 

Returning to (^), let now s'q and s'q be singular points of fg and fa, respectively. Then the integrand in 
has singularities at values s' = s'l and s' = s'2 with imaginary parts given by: 

Im {s'l) — Im {s'q) — e, and Im (s'2) — [Im {s'q) — e] . 

It is then clear that if both s'q and s'q are on or below the real axis, then the corresponding singularities 
of the integrand will be on opposite sides of the real axis iff a/ (a — 1) < 0, i.e., < a < 1. For values of 
a outside this interval, both singularities fall on the same side of the real axis, and the contour of the s' 
integration can be closed with a semicircle that doesn't contain either of them. Thus if all the singularities 
of the functions fQ and fQ occur on or below the real axis, only the interval (0,1) in the integration over 
a contributes to the integral ^ This reasoning is only a slight modification of that presented in Q . The 
integral in then becomes 



^ We are tacitly assuming that the integrals are ultravioletly convergent and thus closing the contour with a semicircle 
at infinity introduces no additional contribution. 
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4^ y_ ^J/'^J T2^""^^ + ^« (-«"'^('?^'' ^)) • (7) 

In dealing with fermion propagators we encounter integrals with powers of the momenta in the numerator 
as shown in (^. As discussed in in such cases the previous reasoning fails for some components of the 
integrals (usually referred to as "bad" components). The issue then arises as to whether or not one can avoid 
all such components. In all such cases we are able to avoid them and compute only the "good" components. 
The "bad" components are recovered by the requirements of Lorentz symmetry. 

We see that to preserve the rules of light-cone theory we must make assumptions about the location of 
the singularities, much like those necessary to justify a Wick rotation. 

A comment on as compared to a typical calculation of a perturbative diagram in a light-cone approach 
is in order. In most light-cone calculations, the integral over (equivalent to the integral over s' in (|^)) 
is absent. An on-shell condition that determines the value of is used, instead of integrating over all 
values. In other words, the calculations are performed in time-ordered perturbation theory (TOPT). This 
has become so customary that the rules of light-cone TOPT are often identified with the rules of light-cone 
theory. The equivalence between TOPT (either light-cone or equal-time) and the Feynman rules is usually 
shown by closing a contour of integration to pick up a singularity in a propagator and thus put a particle on 
shell. This can be done in the equal-time formulation (by doing first the integral over q^ , for example) or in 
light-cone formulation (by doing first the integral over q~). One should refrain from doing so in dealing with 
the Schwinger-Dyson equation. Most certainly, the singularities of the Green's functions involved are more 
complicated than just a pole, for one thing: this is not perturbation theory. While such practice comes in 
naturally in the equal-time formulation, where TOPT is not very customary, it raises a flag in a light-cone 
approach. Hopefully, this comment addresses that flag. 



III. MODELS 



There has been an extensive program of research on the Schwinger-Dyson equation during the past 
decade. See |^,^ for reviews. The usual approach is to model the gluon propagator and to use symmetries 
to express the dressed vertex in terms of the propagators, or to use the free vertex. The model is constrained 
by the vacuum condensates. We follow this general procedure in our light-cone approach. Since our approach 
allows one to define meson form factors for all momentum transfers, which is not true of instant-form 
approaches, applications of our solutions to hadronic properties will be most interesting. 

The Lorentz tensor structure of the gluon propagator has the general form: 

here ^ is the gauge parameter. The choice of a model is to pick a value of ^, such as ^ = for the Landau gauge 
and the form of the function T)(k^). The two essential nonperturbative features that the gluon propagator 
must be consistent with are confinement and the vacuum condensates. It has long been known that with 
the gluon propagator having a behavior there is confinement and one can fit the string tension. 
On the other hand, it is also known that with the instanton liquid model [T5| ] one can fit the quark and 
gluon condensates. The two models that we investigate in the present work are the polynomial model and 
the instanton model, motivated by these two features of nonperturbative QCD. 
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A. Polynomial model 



We work in Landau gauge with ^ = 0, and take as our model gluon propagator 

D,Ak) = - D{k% (9) 



with the parameter x introduced to allow us to use a Feynman-like gauge (the gauge where D^i^ik^ c>c 
Thus X = for Feynman-like gauge and x = 1 for Landau gauge. See Refs [Q,?] for discussions the choice 
of gauge and gauge invariance for Schwinger-Dyson equation in QCD. 

It has been shown elsewhere JlGf that renormalization group arguments yield an approximate relation 
between the renormalized coupling constant, the renormalized gluon propagator and the effective coupling: 

glDRi^) « glffi^) = 47rae//(fc'), (10) 

where the subscript R denotes renormalized quantities, and Qeff is the effective running coupling constant. 
The renormalized coupling is related to the running coupling by g\ = g^jj{Q^) , where is the 

renormalization point. 

As discussed in section 6.1], Eq.(|l^), and the fact that in Landau gauge only the combination g^D (^k^) 
enters the Schwinger-Dyson equation, reduces the problem of modeling the non-perturbative part of the gluon 
propagator to that of modeling the non-perturbative part of cteff- We model ae// with the expression: 

«e//(fc^) = E(-ir^4pT7^) • (11) 

Although the known logarithmic behavior of cteff in the high energy limit cannot be fitted accurately 
with this form, a reasonable approximation can be obtained up to energies well beyond our main region of 
interest. Thus, for example, with = 2, sq = 1 GeV^, Ai = 0.222, ci = 0.07, As = 0.25, C2 = 0.6 the 
values for aeff obtained from our fit are within the accepted error bars as published, for example, by PDG 
Ip^l f\ for energies from about 1.3 GeV to 350 GeV [see Fig |j. Fits good to higher energies can always be 
achieved by adding more terms in (^ij) with smaller q. As can be seen in Fig. ^ with these parameters, our 
model underestimates the accepted values for a^f f in the important energy range of a few hundred MeV to 
about 1.2 GeV, although it is widely believed that it overestimates them in the deep infrared. 
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Numerical data for astrong taken from the PDG at [ittp://www-theory.lbl.gov/ -ianh/alpha/alpha.html 
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FIG. 2. Fitting aeff with Eq.(|ll[). See text for values of parameters, a) For energies up to 350 GeV. b) A close-up 
for lower energies. 



There is little certainty about what happens to the effective coupling or the gluon propagator at very low 
energies. Some results |l8| , |l6| suggest that the pole at fc^ = gets enhanced and could be as strong as 1/fc'*. 
This would provide an explanation for confinement as shown by 'tHooft [Q. Other results ||l^;|22l suggest 
that the pole gets softer or disappears and that the gluon picks up an effective mass. Our model describes 
an infrared-enhanced gluon propagator, although not as strongly as 



B. Instanton model 



As we shall discuss in section VI below, one cannot obtain the known quark condensates with the 
polynomial model. On the other hand, it is known that with a pure instanton model one can obtain the 
quark condensate but not the string tension. For this reason we consider a model with quarks propagating in 
the instanton medium and in addition the gluon propagator having a structure to get confining effects 
of the far infrared. The starting point is the solution for the instanton using the classical action |Q, which 
gives for the instanton color field 

= §^ (12) 



P 

192/ 



(.t2 + p2)4 



where p is the instanton size. The quark zero modes in the instanton background p5[ , for the -I- mode with 
the instanton at position z, are 

= ^ , „ ;-r7r „ 7a (a:^ - Za)U, (13) 

where p is the instanton size and U is a unitary color-spin matrix. From this the widely used model of the 
quark propagating in the instanton-antiinstanton medium pq| was derived. 
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Ai{p) -1 = 

Bi{p) = Kp'fi^p) 

Hp) = ^ - (3/ob) + h{p)) X Xib). (14) 
The quark propagating in the instanton-antiinstanton medium is illustrated in Fig. pi 



FIG. 3. Quark propagating in instanton medium 



In our model with instantons we include this propagator in the Schwinger-Dyson equation by adding it 
to the free propagator. The justification for this model is that the instantons give the nonperturbative 
QCD at the length scale of about 1/3 fm, while the polynomial model gives the far infrared behavior. As a 
simple example of the physics missing if one uses only the polynomial model for the giuon propagator and 
ignores instanton effects consider a delta- function model for the gluon propagator. See, e.g., Ref. [|3[ for the 
solutions to the Schwinger-Dyson equation in such a model. The quark and mixed quark condensates are 
given by 



o r+oo 

g-(0)g(0):|0)=-— dSS 



B{S) 



SA^{S) + B^{S) 

(15) 

■ G(0)g(0) : |0) = 
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r+°o ,oo i n B(S){2-A(S)) iilB{S)[2SA{S)(A{S)-l)+B''(S)] \ 
Jo "'^'^ \^ SA\S) + B-^{S) ^ 16(SA^(S)+B^(5)) /• 

One can fit the strength of the delta function gluon propagator to obtain the correct quark condensate. Then 
from Eqs. (15) it is simple to calculate the mixed quark condensate. The result is an order of magnitude 
different from the phenomenological value. This is an example of the need to include both the long-distance 
and the medium-distance nonperturbative QCD effects. We return to this when we discuss the results of our 
calculations. For completeness we include here the equation for |^, Appendix C] 
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f 



3 
4^ 



dS- 



SM{S) 



A{S)[S + A'P{S)]^ 
Notice that in these equations S denotes Euchdean momentum. 



M{S) 



SdM' 
2~dS 



C. The vertex 



The results reported in this paper were all obtained with the approximation (p, q) = 7^ (rainbow 
approximation). The popularity of this approximation is justified by its simplicity. It suffers the serious 
drawback of violating the Ward-Takahashi identity (WTI): 

{p-qrr^ip,q) = S-\p)-S-\q) (16) 

This identity holds in QED, and the forms for the vertex discussed below were proposed for QED. Such forms 
are, however, often used for QCD as well, since the analogous identity for QCD (known as the Slavnov- Taylor 
identity) reduces to (jl^) if ghost effects are ignored. Ignoring these effects is believed to be important in the 
infrared region. 

This identity can be used to express the longitudinal part of the vertex in terms of the propagator. Ball 
and Chiu |2ll have proposed the form: 



r,{p, q) = r^^(p, q)+Y.r (p^ • q) T;(p, q), (17) 
j-1 



where the T^'s are eight transverse tensors (i.e. tensors that satisfy {p — q)'^T^{p, q) = 0, and thus don't 
contribute to the WTI). The form of these tensors is given in Eq.(3.4) of [^. The eight scalar functions 
are not constrained by the WTI. T^'~^ is given by: 

r^l^^p) = 7. {A{q')il /3i) + A(/)/3i) + 

(18) 

+(. + P). - P.)) ^(4^ - (.+p) /(4-y . 

_ pz q — p 

This is really a small modification of the vertex proposed by Ball and Chiu. In their form /3i = 1/2, but 
other values are often used to explore the effects of the transverse part. Besides compromising the symmetry 
of , changing the value of (3i can have serious effects on renormalizability, namely because of its effect on 
the transverse part, as discussed below. The WTI itself is satisfied for any value of /3i, and for any choice of 
/"s. The converse is also true: any vertex that satisfies the WTI can be expressed through (|l7| , [l8| ). 

The vertex has been further constrained by other studies, in particular. Burdens and Roberts [ p2[ list 
a number of requirements should satisfy, besides the WTI. See also |Q, section 3.7]. Exploiting these 
requirements, and particularly the need to maintain multiplicative renormalizability in QED, Curtis and 



Pennington ||23| narrowed down the last term in Eq.(17) to 



(.,.)= /3.(7.(. -P)- i,+p)M- ^)) ^ (M^(,^) +iiP(p^))- - ^''^ 
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Here too /?2 = 1/2 in the original formulation. Applying the approach explained in section^ we have found 
that when using 



one must set consistently Pi — [32- This ensures, for example that the function A remains free of divergencies 
in Landau gauge, a well established result. Notice that this implies that using V^'-^ by itself is, in general, 
inconsistent. 

We are currently investigating using this form in our model. Although that entails neglecting the effect of 
ghost fields in the STI, it is a significant improvement over the rainbow approximation. 



IV. INTEGRATING THE EQUATIONS. 

In order to solve (Q) in our model, we need to calculate integrals of the type (here a/ = 1 + q): 




(20) 



{ Coo (P^ /q , aO , P'' Cio (p2 ; /q , aO , P^p"" C20 (p^ /q , a; ) + 5^"^?oi (p' ; /q , a; ) } , 



where we have exploited the Lorentz structure to define the scalar quantities C, 
These integrals include as a particular case the well known integrals: 




(21) 



{/oo(p^; ao, ai),P^Iw{p^] oo, ai),p^'p'' hoijP'; ao,ai) + g'"" Ioi{p'^; oq, a;)} . 



For (r,j) = (0,0), (1,0), (2,0), (0,1), the I's are known to be: 




X 



(22) 




(a(l - a)(-.s) + (1 - a)ml + amj - ie)'^o+ai-d/2-j ' 



^'^'+'^-1(1 - a} 



where we have used the compact notation 



r 



01,02, . 
61, &2, • 



, a„ 1 r(ai)r(a2) . . .r(a„) 

J " r(6i)r(62)...r(&„) ' 



For the C's, we proceed as explained in section O to get: 
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Crj{s, fQ,ai) = r 



da 



a; + 1 - rf/2 - j 
ai 



+00 



(23) 



(a(l - a)(-.s) + (1 - a)s' + Q;(mf - if)Yi+i-d/2~j ■ 



then (H) agrees with (pi). This 



At this point we must make sure that if we set /q (s') = 1/ (s' — "^o)""; 
is most easily done by closing the contour of integration of the s' variable with an infinite semicircle in 
the lower half of the complex plane to pick up the singularity in fq. This procedure, however, uses the 
analytic form of fq. We need a procedure that would rely only on the numerical values of fq. We perform 
the integrations by closing our contour with an infinite semicircle on the upper half of the complex plane to 
pick up the singularity (in s') of the expression on the second line of (p3|). For mi = 0, the I' s reduce to ^ 

Irj{s] ao, ai) = 



j(_l)-("o+a,)(_l/2p 

(47r)''/2(mg)'^o+'^i-''/2-J " 



oo + a; — d/2 — j, ai + r, d/2 + j — ai 
ao,ai,d/2+j + r 



For the C 



+r 



1 + 6 - a 

b, c — b 

c, 1 — a 



2^1(00 + ai~ d/2 ~ j, ai + r;d/2 + j + r; s/ml). 
(-s)-* / ds7Q(s')(-s')''""2Fi(l + 6 - c, 6; 1 + 6 - a; s' / s) + 



ds'fQ{s'){-s') ''2Fi{a,b-c;s/s') 



s <0 



(24) 



b, c — b 

c, 1 — a 



s = 



(25) 



c- 6 
1 + c — a — 6 



+r 



b, c~ b 

c, 1 — a 



- / ds7Q(s')(s - s'Y-''^\s'f-\F^(l - a, 1 - b; 1 + c - a - &; 1 - s/s') + 
f^s'/Q(s')(-s')~"2-Fi(a, fe; c; s/s ) s > 0, 



''if ao is not an integer, then this singularity will be a branch point with its corresponding branch cut, as opposed 
to a pole. 

^The expansion of this expression about d = 4 is probably more familiar to the expert in dimensional regularization. 
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where a — ai + 1 — d/2 — j, b — ai + r, c — d/2 + j + r. 



V. SAMPLE CALCULATIONS IN MINKOWSKI SPACE. 



In this section we discuss sample calculations in Minkowski space. The ability to work in Minkowski 
space is one of the important features of our light-cone formulation of the SD equation, in this section, 
however, we only use polynomial models, which allow convenient Minkowski space forms. As we shall see, 
the infra-red properties in these models are not consistent with the known condensates, and solutions with 
instanton contributions to the gluon propagator are much more satisfactory. Since the instanton forms have 
been obtained in Euclidean space, however, there is no unambiguous Minkowski space formulation. We 
discuss our results with models including instantons in the following section. 

Eq.(p5|) is a generalization of the well known result Eq.(24) to a wider class of functions /q. The light- 



cone approach used, and in particular the treatment it makes of the singularities, seldom goes without some 
aftertaste. Here is a well known example where using light-cone variables leads to a dead end ^ 

f dkodki 
lim / —Tj -75 ^ — n > 1 

e^a+ J {kl -kl^ml^ le)" 

This integral is perfectly convergent and can be computed, for example, by closing the contour in the fcp- 
integration with a semicircle at infinity. The integral is, of course, non- vanishing. Using light-cone variables, 
one would have the integral: 

f dk+dk- 

/ „/, , , ^ n > 1 

e~>o+ 7 2{k+k- -ml + ieY 

If we now think of closing the contour of, say the /c^-intcgration, we see that the singularity occurs at 
(ttiq — ie)/k~^. But then the singularity can always be avoided for any k~^ ^ 0. Only for fc+ = is the 
^"-integration non- vanishing, but then, it must have a delta function type singularity. 



The purpose of this section is test Eq.(25) against such anomalies, as well as to point out some important 



technicalities. We input a "first guess" for the quark propagator: fojs') = l/(s' — ttiq) with ttiq = 2 GeV, 



and use the model for the gluon propagator described in section [II 



k'^ + ie J \k^ + ie 

where sq = 1 GeV^, ci = 0.07, Ai = 0.222, C2 = 0.6, A2 — 0.25. We also discuss the ultravioletly divergent 
case ci = 0. For the vertex we set =7^. We use Eq ( p5| ) to compute the functions A and B with this 
input. We exploit the fact that for the first iteration, the integrals could also be computed analytically 
using Eq.(^J). We present graphs comparing the two outputs. We stress that these are not solutions to the 
Schwinger-Dyson equation but merely the output of running the iterative procedure once. Solutions to the 



Schwinger-Dyson equation are presented in section VI 



O.L. wishes to thank Dr. Matthias Burkardt for bringing this to his attention. 
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FIG. 4. The quark propagator after one iteration compared with the analytic result, a) analytic result (|24|); 
numerical integrations of (E5|):6^ with tail contribution, c) without tail contribution. Notice the appearance of a 
singularity and an imaginary part. 



Keeping with tradition, the graphs are shown with a reversed x-axis. We use the variable Se = — s, 
equivalent to the Euclidean invariant momentum squared. We show, however, negative values of Se^ which 
lay outside the usual Euclidean space, and represent the time-like region. The numerical integration is 
carried out up to Semax = 64 GeV^ and an estimate of the contribution of the remaining integration out 
to infinity ("tail contribution") is added. As will be seen on the graphs the consequences of omitting this 
tail contribution are numerically important. 

In fact, in order to reproduce the numerical values of Eq.(]2^), Eq. (p5|) demands an elaborate integration 
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procedure. On one hand there is the ever present fact of having to integrate over all energies. Renormal- 
izability tells us that there is much freedom in dealing with the very high energy region (it can even be 
dropped, as in cut-off renormalization schemes), and whatever numerical effect that has on the integrals can 
be compensated by readjusting a finite number of parameters. From this point of view the effect of omit- 
ting the tail contribution, although numerically important, should be of no physical consequence. Eq.(24), 



however, represents the most accepted treatment of the very high energy regime: the integrals are carried 
out to infinity without the introduction of cut-offs or form factors that could spoil the symmetries of the 
theory. Divergencies, when they exist, are treated with the dimensional regularization method. In order 
to have Eq.(|24|) as a reference, we choose to carry out the integral in Eq.(p5|) out to infinity as well. In a 
numerical procedure this means adding an estimate of the tail contribution. On the other hand, we have 
a "floating singularity" in the integrand (a singularity at Se = S'^, as opposed to one at a fixed point on 
the grid). This is particularly important for the infrared-enhanced terms (those with larger exponent q). 
Given the strength of the singularity at Se = S'^, the contribution from the region Se ~ S'^ is important. 
Sampling enough values of S'^ immediately to the left of Se, is impossible when Se ~ Semin^ unless the 
grid is refined each time (the analogous problem around Semax g^ts handled by the addition of the tail 
contribution). This can contribute to inaccuracies in the very low energy region, a most undesirable effect. 
Refining the grid with each iteration, however, can easily lead to exponential growth of the computational 
time. 
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FIG. 5. The numerically renormalized quark propagator after one iteration in an ultraviolet divergent case, com- 
pared with the analytic MS result, a) analytic resul t (p^ ; numerical integrations of (p5[): b) with tail contribution 
and MS subtraction, c) without tail contribution or MS subtraction. 



One solution is to use a model propagator with no infrared enhancement, and that falls off rapidly at large 
energies, e.g., a Gaussian model. We don't feel, however that Gaussian models are appropriate in the time- 
like region, and one of the purposes of this paper is to explore the effects of a tunable infrared enhancement. 
Our solution is to sample numerically only the part of the integrand involving the fermion propagator and 
the vertex, and to treat the part coming from the boson propagator (and containing the floating singularity) 
analytically. In the language of numerical methods, it is much along the lines of Gaussian quadratures. With 
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this, and the addition of the tail contribution, the procedure becomes accurate, and reasonably fast. 

The graphs in fig. |^ show the real and imaginary parts of the functions A and B after one iteration using 
the parameters above. As can be seen from the graphs, the procedure is quite accurate (provided the estimate 
of the tail contribution is added!). The relative error is of the order of thousandths of a percent, although 
it is somewhat larger right before (but not after!) the singularity at Se — —niQ. The imaginary parts are 
obtained by simply applying the correct ie— procedure around the singularity in /q. 

Figure |^ shows the renormalized real parts of the same functions, with the crucial difference of setting 
ci — 0, rather than ci = 0.07. This causes an ultraviolet divergence. The imaginary parts in this calculation 
contain no divergence, so they remain basically the same as in the convergent case, and thus are not shown. In 
this case, the analytic results are renormalized in the MS scheme. The numerical calculation was performed 
with d ~ 4: + 2e, where e = —10"^. The divergence comes in the tail contribution (which is almost constant 
for the divergent term), and it is canceled by the 1/e subtraction used in the MS scheme. The well known 
result that A is free of ultraviolet divergences in Landau gauge renders this function insensitive to whether 
or not the subtraction is performed. Addition of the tail contribution still significantly affects the non- 
divergent term. Even for B, it is possible to obtain a finite result by just dropping the tail contribution 
(cut-off renormalization) , but we choose to integrate to infinity and perform the 1/e subtraction. We do not 
include in this paper solutions to the Schwinger-Dyson equation involving this renormalization procedure, 
but thought it appropriate to show how such solutions could be obtained. We point out that dimensional 
regularization is not believed to give physically correct results in non-perturbative calculations. This method 
could be used to explore that question. 



VI. RESULTS WITH INSTANTONS CONTRIBUTIONS 



In the present section we treat models that contain instanton contributions. Since there are a variety 
of possible models, we present results obtained with different approaches. On one hand, we want to study 
Eq. (|l]) in its original form. On the other hand we want include the effects of the propagation in the instanton- 
antiinstanton medium, contained in (|l^). We compare three different alternatives for including the effects 
of (|l^) in (|l|). All in all we study four different equations: 

S-\p)^S^'(p)-j:gj^{S{p)) (26) 



S-\p)=Sy\p)~J:sj^{S{p)) (27) 
S~\p) = So\p) -J:sD i^(P^ + ^Hp)) (28) 

S-\p) - ST\p) -^SD (^iP^ + ^^Hp)) , (29) 

where 5*0 is the bare propagator, Sj is the propagator given by (|lj), S is the full dressed propagator, and 
'^SD is the quark self-energy, thought of here as of a function of S{p). The subscript SD is introduced here 
to distinguish this self-energy from the kernel of a different equation, discussed below. Thus, Eq. (p6|) is just 
Eq.(|l]) and does not explicitly include the instanton effects contained in Eqs. (p^), while Eqs .(^7| , |2^, |29| ) 
propose three different ways in which these effects could be included. The form of the function "^JgjjXSip)) 
(see Fig. |^ and Eq.(|l|)) depends on the model chosen for the gluon propagator, the gauge and the vertex. We 
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use two different sets of parameters for tfie gluon propagator and the rainbow approximation for the vertex. 
All calculations were performed in Landau gauge and in the chiral limit {rric = 0). The parameters for the 
gluon propagator are: 



set a. I ^ Q 25 C2 = 0.6 ' ™^ ^ \ As = 1.5 C2 = 0.85 

In the table below we compare the different solutions in terms of the values they yield for the quark 
condensate, the mixed quark condensate, /^r, and the value of AI = B/A at = Fl. 



Vertex: — j'^ Rainbow approximation 
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"see footnote |^ at the bottom of page 



Most of the results obtained with Eq. (|2^) (no instanton effects) set a, appear too low. This is most likely 
due to the fact that the corresponding coupling is too weak in the crucial region of a few hundred MeV to 
about 1.2 GeV. Inclusion of the instanton effects seem to improve matters considerably, perhaps with the 
exception of Eq. (psl) . 



A plausible explanation for this maybe as follows: the solution (14) was also obtained from a self consistent 
equation |26|, which we write symbolically as: 



One could think of including both types of one-particle-irreducible diagrams into the self-energy: 



(30) 



(31) 



^This value is often compared to the constituent quark mass (Mq) of the non-relativistic quark model (NRQM). 
As we show in Figs. |^ and ^, the solutions to the Schwinger-Dyson equation vary quite rapidly upon entering the 
time-like region. The value Mq should therefore, at best, be regarded as a rough guide to the value of B/A at = 0. 
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Solving ( pO| ) for Sq ^ and substituting into ( |3l| ) 

= sy\p) - (^b)) - Ej (^(P) - ^/(P)) (32) 



2 
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FIG. 6. Two different fits to a, 



Eq. ( p2D contains the effects of the quark propagating in the instanton medium as suggested in section [II . 
An important feature is the presence of 5"^^^ instead of 5*0"^ in the inhomogeneous term, which occurs also 
in Eqs. ( p7| , p9| ), but not in Eq. (|2^) . In fact, Eq.(p7|) is an approximation to Eq.(^ up to terms of order 
Y^j {S{p) — Si{p)), which can be considered small, as opposed to Sq^ — SJ^, which obviously is not small. 
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FIG. 7. Solutions to Eq.(pq) using set a. The left part of the graphs {Se < 0) represents the time-like region. 



The results obtained with Eq.(p6|) (no instanton effects) set b, appear closer to the correct values, often 
overestimating them. This is most likely due to the fact that set b contains stronger infrared enhancement 
and overestimates the coupling the intermediate region. A comparison of both sets in terms of how they fit 
the accepted values of a^ff is given in Fig. ^. Inclusion of the instanton effects with set b, in most cases, do 
not lead to worse overestimation, as could be naively expected. With set a we have a situation where the 
polynomial part of our model is contributing very little in the intermediate region (which we believe to be 
crucial to the quantities we are computing), and the instanton effects bring in the bulk of the contribution. 
With set b both parts are contributing significantly, but their effects do not add. Thinking along the lines of 
Eq.(|l]), this probably means that the two types of diagrams are of different nature and not always interfere 
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constructively. This suggests that the instanton effects are not being double-counted with this procedure, 
although such conjecture needs further testing. 
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FIG. 8. Solutions to Eq.(E6|) using set b. The left part of the graphs {Se < 0) represents the time- like region. 



Figures ^ and ^ show graphs of the solutions obtained using set a and set b, respectively, and Eq. ([2^). 
The focus is on the low energy space-like region (right half of the graphs) and the cross-over to the time-like 
region (left half of the graphs). Already at small time- like energies (the graphs show energies up to about 
360 MeV for set a and 1.3 GeV for set b) the behavior of the functions has changed dramatically from their 
behavior in the space-like region. The function A, for example, stops growing; and, say, the function M for 
set a, (see Fig|7|) which has barely risen from its asymptotic value of M{p'^) rUc = at large space-like 
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to about 90 MeV at = 0, is already close to 400 MeV at w {mQMeVf. 
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FIG. 9. Solutions to the Schwinger-Dyson equation using set a in the space-like region, a) Result Eq.(14) from 
Ref[26], b) Eq.(26), c) Eq.(27), d) Eq.(28), (e) Eq.(29). In (a) A = 1. 
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FIG. 10. Solutions to the Schwinger-Dyson equation using set b in the space- like region, a) Result Eq.(14) from 
Ref[26], b) Eq.(26), c) Eq.(27), d) Eq.(28), (e) Eq.(29). In (a) A = 1. 



Of interest is also the behavior of the function D (Se) = SeA"^ (Se) + {^e), the denominator of the 
propagator. The vanishing of this function would indicate a pole in the quark propagator. Extrapolation of 
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the behavior in the space-hke region would suggest a pole at much smaller time-like energies than those we 
have plotted, particularly for set a (for set b some turning around is already noticeable at very low space-like 
energies). The behavior of D right around its zero, if there is one, would be difficult to extract, since the 
functions A and B become singular there. In fact, D seems to drop fast just to the left of the regions we have 
plotted. Thus, there seems to be a pole in the quark propagator on the real axis, in our polynomial model. 
The pole would be at « (360 MeV)^ for set a, and ~ (1.3 GeV)^ for set b. Further investigation of the 
solutions deeper into the time like region, is necessary before this can be ascertained. Notice that this refers 
to Eq.(|^ only. 

Figures ^ and |l^ compare the solutions to Eqs.(|27|, [28| [2S| ) for set a and set b, respectively. They include 
only the space- like region (Euclidean space), since they use the result (|lj) from which was obtained in 
Euclidean space. The analytic continuation of result (^^ into the time-like region exhibits a branch point 
at = and thus has not been used. 



VII. CONCLUSIONS. 



We have shown how to obtain a light-cone form of the Schwinger-Dyson equation for the quark propa- 
gator. For QCD the key question is the infra-red behavior of the gluon propagator, for which we have used 
two models. Knowing that the 'tHooft model gives confinement, we have used polynomial models for the 
gluon propagator based on that model. Although the instanton liquid form for the gluon propagator does 
not confine, and therefore does not have the correct far infra-red behavior, it provides the main mid-range 
QCD interaction, and allows fits to the condensates, which is essential for obtaining hadronic properties. 
With our approach we can obtain solutions for both the enhanced infra-red behavior of the polynomial type 
models or the regular infra-red behavior of the instanton model and recent lattice gauge calculations [p2| . 
With the polynomial models we are able to work in Minkowski space and obtain solutions. With the models 
that include instanton effects we are able to obtain solutions that give much better agreement with the 
phenomenological values of the condensates. 

We consider the present work is exploratory. It provides the framework for obtaining light-cone QCD 
propagators that can be used to obtain light-cone models of hadronic BS amplitudes for studies of hadronic 
properties at all momentum transfers. 
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